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ABSTRACT 

The relatively nearby (distance=24.1 Mpc) elliptical galaxy NGC821, hosting 
a central massive black hole but inactive at all wavelengths, was observed with 
Chandra for a total exposure of 230 ksec, to search for nuclear emission and gas 
available for accretion. Within its optical image, 41 sources were detected, with 
spectral properties typical of low mass X-ray binaries (LMXBs). The fractions 
of LMXBs in the field and in globular clusters were determined, together with 
their X-ray luminosity function (XLF) down to L(0.3-8 keV)=2xl0 37 erg s _1 . 
At the galactic center a source of L(0.3-8 keV)=6 x 10 38 erg s _1 was detected for 
the first time, slightly extended. Its spectral shape is quite hard (T = 1. 49^13), 
without intrinsic absorption. It is surrounded by three sources with spectral 
shape typical of LMXBs and luminosities on the brightest end of the XLF. One 
is consistent with being pointlike; the others could be the superposition of few 
point sources and/or truly diffuse emission, with one resembling a jet-like feature. 
Diffuse emission was detected out to R ~ 30", and comes mostly from unresolved 
LMXBs, with a minor contribution from other types of stellar sources. Different 
lines of investigation consistently provide no evidence for hot gas. Hydrodynami- 
cal simulations show that stellar mass losses are driven out of NGC821 in a wind 
sustained by type la supernovae, but also hot accreting gas within a very small 
inner region. A companion paper presents further observational results from 
Spitzer and the VLA, and possible accretion modalities for this central massive 
black hole. 
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1. Introduction 

NGC 821 is an isolated elliptical galaxy at a distance of 24.1 Mpc (Smith et al. 2004; 
see Tabled]), with very regular optical isophotes of disky shape (Bender et al. 1994, Lauer 
et al. 2005) and an old and metal rich stellar population as typical of elliptical galaxies 
(Proctor et al. 2005). No cold (HI) or dusty ISM has been revealed in it (Lauer et al. 2005, 
Sarzi et al. 2006). NGC 821 is also one of those ~ 30 nearby galactic spheroids for which 
a central supermassive black hole (MBH) has been claimed based on resolved dynamical 
studies (Ferrarese & Ford 2005). The mass of this MBH is 8.5 x 1O 7 M (Table ED and its 
Eddington luminosity is L^dd ~ 1-1 x 10 46 erg s _1 . However, this MBH is extremely quiescent 
at all wavelengths (see also Pellegrini et al. 2007), and represents an excellent example of the 
class of "inactive" nuclei, that are very common in the local universe (e.g., Pellegrini 2005a, 
Ho 2005) and key to our understanding of the mechanisms of AGN evolution and nuclear 
feedback, that have been fundamental for shaping our universe (e.g., Springel, Di Matteo 
& Hernquist 2005). Since it is nearby, NGC821 is a prime target for observing a quiescent 
nucleus and its surroundings. 

X-ray emission is a key symptom of nuclear activity resulting from accretion onto a 
MBH (Rees 1984), and the hot ISM that could provide a source of fuel is readily visible 
in the X-rays (e.g., Fabian & Canizares 1988; Loewenstein et al. 2001; Pellegrini 2005a). 
Previous observations of NGC821 with the ROSAT PSPC had placed an upper limit of 
2.8 x 10 40 erg s^ 1 (O'Sullivan et al. 2001) on the total X-ray luminosity, suggesting a low 
content of hot gas (e.g., Kim et al. 1992). With the sub-arcsecond resolution of Chandra 
we could attempt to resolve the different components of the X-ray emission, i.e., the hot 
gas, the nucleus and the population of low mass X-ray binaries (LMXBs; Fabbiano et al. 
2004, hereafter F04). This previous shallow (39 ks) Chandra observation, performed in 
2002, revealed diffuse emission in the central galactic region, possibly from a hot interstellar 
medium. We did not detect any source that could be unequivocally identified as the galactic 
nucleus (for which we derived L x /L Edd < 10~ 7 ), but an intriguing S-shaped feature crossing 
the center of the galaxy was observed. Explanations for this feature included either a weak 
two-sided X-ray nuclear jet, or a hot gas filament, at a temperature higher than that of the 
surrounding gas. This opened the possibilities that the accretion power could end mostly 
into mechanical rather than radiative power (Di Matteo et al. 2003, Pellegrini et al. 2003a, 
Fabbiano et al. 2003), or that the accretion flow may be disrupted by nuclear feedback (e.g., 
Ciotti & Ostriker 2001, Omma et al. 2004; see also Soria et al. 2006a,b). 

In this paper we report the results of the analysis of deep Chandra observations (for 
a total exposure of nearly 230 ks) aimed at defining the properties of the various sources 
in the central galactic region. With the deep Chandra exposure, we can constrain all the 
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components (nucleus, hot gas, unresolved sources) contributing to the emission, a necessary 
and critical step when the goal is to understand the accretion process. For example, we can 
detect a significant portion of the LMXB population, which was barely visible in the 2002 
data, and thanks to the calculation of its X-ray Luminosity Function (XLF) we can now 
discriminate between the contributions to the diffuse emission of the unresolved LMXBs and 
of the hot ISM (e.g., Kim & Fabbiano 2004). Moreover, determining the properties of the hot 
ISM also gives clues on the way the ISM evolves in this galaxy, and more in general in galaxies 
of similar optical luminosity (e.g., Pellegrini & Ciotti 1998, Sansom et al. 2006, David et al. 
2006). A companion paper (Pellegrini et al. 2007) presents an observational campaign aimed 
at determining the properties of the nuclear emission of NGC 821 at various wavelengths, 
including more sensitive proprietary VLA and Spitzer IRAC observations, together with 
archival HST observations. 

The paper is organized as follows: in Section [2] we describe the Chandra observations 
and data preparation; in Section [3] we present the X-ray properties of the whole LMXB 
population of the galaxy, also making use of HST WFPC2 images to identify those in field 
and in globular clusters; in Sections HO and [5] we report the results of the X-ray analysis 
respectively for the nuclear region and the diffuse emission; in Section [6] we discuss the 
implications of our observations; in Section [7] we summarize our results. 

2. Chandra observations and data preparation 

NGC 821 was observed with Chandra ACIS-S (Weisskopf et al. 2000) seven times 
between November 2002 and June 2005, for a total exposure of 230 ks (Table [2]). Here we 
consider both the archival observations (ObsID 4408 and 4006 in Tab. [2]), already presented 
in F04, and the new unpublished data obtained in 2004 and 2005. 

The satellite telemetry is processed at the Chandra X-ray Center (CXC) with the 
Standard Data Processing (SDP) pipelines, to correct for the motion of the satellite and 
to apply instrument calibration. The data used in this work were reprocessed in custom 
mode with the version 7.6.0 of the SDP, to take advantage of improvements in processing 
software and calibration, not available at the time of the original processing. Verification 
of the data products showed no anomalies. The relative astrometry of the different obser- 
vations was corrected by using detected sources (above a 3cr threshold) as reference points 
( [http://cxc.harvard.edu/cal/ASPECT/align_evt7l . 

The data products were then analyzed with the CXC CIAO v3.0.1 software and the 
HEASARC XSPEC package. CIAO Data Model tools were used for data processing, such 
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as screening out bad pixels and producing images in given energy bands. Calibration files in 
CALDB version 3.1.0 of June 2005 were used. No significant background flares were observed 
in these data, so no further screening was necessary. A time-dependent gain correction 
(Vikhlinin et al. 2003) was applied to the SDP Level 2 event files of the first two observations 
before further analysis. This correction was already included in the SDP Level 2 event files 
of the other observations supplied by the CXC. The seven observations were then coadded 
to produce a deep image, using the CIAO task merge jail, which takes as input the event 
files and the relative aspect solution (ASOL) files, and reprojects the coordinates of all the 
observations onto the first one. 

Below we report the results of the analysis of these data. 

3. The Low Mass X-ray Binaries population 

Figure Q] shows the merged X-ray image, with the D25 ellipse of NGC821 overplotted 
(this ellipse is the isophote of the 25.0 B-mag/square arcsec brightness level, Tab. [1]); also 
marked are the X-ray sources detected by the CIAO task wavdetect. Of the 104 sources 
detected within the ACIS-S3 CCD, 41 lie inside the D25 ellipse. Except for minor contam- 
ination from interlopers (see Sect. 13.21 and 15. 6p . and four particular sources in the central 
region (discussed in Sect. H] below), these sources are pointlike and represent the LMXBs 
population of NGC821 (see Fabbiano 2006 for a review on X-ray binary populations of 
galaxies) . 

For the pointlike sources within D25, counts were extracted within a circle centered on 
the wavcfetect-determined source position, and background counts were estimated locally in 
an annulus surrounding the source. We chose the source extraction radius to be the 95% 
encircled energy radius at 1.5 keV (that is varying as a function of the off-axis angle), with a 
minimum of 3" near the aim point. Similarly, the background was estimated for each source 
from an annulus surrounding it, with inner and outer radii of 2 and 5 times the source 
radius respectively. When nearby sources are found within the background region, they 
are excluded before measuring the background counts. Net count rates were then calculated 
with the effective exposure (including vignetting) for both the source and background regions. 
Errors on counts were derived following Gehrels (1986). When the source extraction regions 
of nearby sources overlap, to avoid an overestimate of their source count rates, we calculated 
the source counts from a pie-sector, excluding the nearby source region, and then rescaled 
them based on the area ratio of the chosen pie to the full disk. Once the correction factor 
is determined, the same factor can be applied to correct counts in all energy bands. For a 
small number of sources which overlap with nearby sources in a more complex way (e.g., 
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overlap with more than 2 sources), instead of correcting the aperture photometry, we used 
the source cell determined by wavdetect to extract the source counts. 

The range of net detected counts for the pointlike sources is ~ 8 — 300, corresponding to 
0.3-8 keV luminosities of 2 x 10 37 - 8 x 10 38 erg s _1 , when using the energy conversion factor 
described in Sect. l3.2l below. Most of the sources are too faint for a detailed spectral analysis, 
therefore their hardness ratio and X-ray colors were calculated in order to characterize their 
spectral properties. The X-ray hardness ratio is defined as HR = (H— S) / (H+S), where S 
and H are the net counts in the 0.5 — 2.0 keV and 2.0 — 8.0 keV bands respectively. Following 
the prescription of Kim et al. (2004), the X-ray colors are defined as C21 = log(Cl/C2) and 
C32 = log(C2/C3), where CI, C2 and C3 are the net counts respectively in the energy bands 
of 0.3-0.9 keV, 0.9-2.5 keV and 2.5-8 keV. These counts were corrected for the temporal 
QE variation, referring them all to the first observing epoch (Nov. 2002, Tab. EJ), and for 
the effect of the Galactic absorption (using the value in Table [I]). By definition, as the 
X-ray spectra become harder, the HR increases and the X-ray colors decrease. For faint 
sources with a small number of counts, the formal calculation of the HR and colors often 
results in unreliable errors, because of negative net counts in one band and an asymmetric 
Poisson distribution. Therefore, we applied a Bayesian approach developed by Park et al. 
(2006), which models the detected counts as a non-homogeneous Poisson process, to derive 
the uncertainties associated with the HR and colors. 

The main properties of the point sources detected in the merged observation within the 
D25 ellipse are summarized in Table [31 where the sources are listed in order of increasing 
distance from the galactic center. Column (1) gives the IAU name (following the convention 
that the name should be "CXOU Jhhmmss.s+/— ddmmss"); columns (2) and (3) give the 
J2000.0 source position and column (4) its associated uncertainty at the 95% confidence level, 
determined based on the source net counts and off-axis angle, from the empirical formula 
given by Kim et al. (2007); column (5) gives the projected distance R from the galactic 
center; columns (6) and (7) give the net counts and the la error; columns (8) and (9) give 
the flux and luminosity in the 0.3-8 keV band, derived using the energy conversion factor 
described in Sect. 13.21 below; columns (10)-(15) give the hardness ratio HR, the C21 and 
C32 colors defined above, together with the respective error determined at a la significance 
level. 



3.1. Spectral properties and optical identification 

Figure [2] shows the C32 - C21 distribution of the point source population; the grid 
gives the location of spectra described by power laws with given values of the photon index 
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T and intrinsic column density Njj. In agreement with previous spectral studies of LMXB 
populations (e.g., Kim, Fabbiano & Trinchieri 1992; Irwin, Athey & Bregman 2003, Fabbiano 
2006), most well-defined colors fall near a typical T = 1.5 — 2.0 spectrum, with no intrinsic 
absorption. 

Three of the sources in Tab.[3]have been detected with a number of net counts sufficiently 
high for a specific spectral analysis (they are marked in Tab. |3]). We performed such analysis, 
extracting from the merged event file the spectra for these sources, using an extraction 
radius of 2"; this radius was chosen to avoid contamination from nearby fainter sources. 
Instrumental response files, weighted by the exposure time of each individual observation, 
were also generated; the background spectrum was derived from a source-free region outside 
the optical body of the galaxy. The results of the spectral analysis (Tab. HI Fig. [3]) confirm 
that these sources are well described by power laws with low-to-moderate intrinsic absorption, 
as typical for the LMXB population. Their intrinsic 0.3-8 keV luminosities are 4.1, 6.9 and 
8.9 xlO 38 erg s -1 . 

With the aid of HST WFPC2 images, we then investigated possible associations of our 
detected point sources with LMXBs belonging to the galactic field or to globular clusters 
(GC). Kundu & Whitmore (2001) studied the globular cluster population of NGC 821 using 
HST WFPC2 images in the V (F555W) and / (F814W) filters; these optical images cover a 
large fraction (~ 60%) of the D25 ellipse of NGC821. Kundu & Whitmore (2001) identified 
105 GC candidates, whose positions, magnitudes and colors were kindly provided to us by A. 
Kundu (2007, private communication). This list of GCs was cross-correlated with the X-ray 
source list of Tab. [31 in order to look for GC/X-ray source coincidences. After correction 
for a small shift of the optical image (see Sect. 14.21 for more details), 6 matches were found 
well within the 95% uncertainty radius of the X-ray position; another potential match has 
quite a distorted shape and uncertain magnitude in the optical (CXOU J020817. 7+105907), 
since it is located at a chip edge, and another one has an offset of 1.24 times the X-ray 
positional uncertainty. These 8 sources are marked in Tab. [3] (see also Fig. [8]). From a 
V-(V — /) diagram of the whole sample of 105 GC candidates, one can see that LMXBs are 
preferentially found in the brightest clusters (4 matches reside in the 5 brightest GCs, and 
6 in the 12 brightest ones), following a general trend already reported (e.g., Sarazin et al. 
2003, Kim E. et al. 2006) that is explained as the probability of hosting a LMXB increasing 
with the GC luminosity. The LMXBs of NGC821 also tend to be associated with the redder 
GCs, in agreement with previous works (e.g., Kundu et al. 2007), even though the small 
number of sources here does not make this result statistically significant. 

From the HST WFPC2 images above, in the overlap region with the D25 ellipse, 17 
X-ray sources have no optical counterpart, at a 50% completeness limit of 24.1 V-mags (A. 
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Kundu 2007, private communicatioro); these can be considered field LMXB candidates, and 
are marked in Tab. [3j In the following we further limit the analysis to sources located at 
radii larger than 10", to avoid possible completeness problems near the galactic center (see 
Sect. 13.21 below). Then, it turns out that ~ 30% of the X-ray sources lie in GCs. From 
Chandra studies this fraction is known to increase with the GC specific frequency of the 
galaxy, or the morphological type (e.g., Sarazin et al. 2003), and goes from ~ 20% in SOs 
to 30-50% in Es (e.g., Fabbiano 2006 and references therein). The fraction estimated here 
for NGC821 is therefore close to what expected, given also its morphological type (E6). The 
fraction of GCs that host an X-ray source is ~ 6.6%, a value broadly consistent with what 
typically observed previously from Chandra studies (3-5%, Sarazin et al. 2003, Kim, E. et 
al. 2006; see also Kundu et al. 2007), but somewhat higher, likely because of the larger depth 
of our X-ray observation compared to the previous analyses. Finally, there is no statistically 
significant difference between GC and field LMXBs. The present sample of X-ray sources 
is too small to investigate other statistical properties, such as the X-ray spectral difference 
between LMXBs in red and blue GCs or the role of galactocentric distance. 



For the point sources falling within the D25 ellipse we estimated an XLF, that in dif- 
ferential form can be expressed as: 



To construct the XLF, we calculated fluxes and luminosities of the point sources in the 
0.3-8 keV band with an energy conversion factor (ECF) corresponding to an assumed power 
law spectral shape with Y = 1.7 and Galactic N H (see Fig. [2] for a justification of this 
assumption). The ECF was calculated with the arf (auxiliary response file) and the rmf 
(redistribution matrix file) generated for each source in each observation. As done for the 
colors, we took into account the temporal QE variation by calculating the ECF in each 
observation and then taking an exposure-weighted mean ECF. The ECF over the 0.3-8 keV 
band varied by ~ 2.5% between 2002 and 2005. We then corrected the XLF to eliminate 
biases and incompleteness effects, which could affect the lowest luminosity range, causing an 
artificial break (Kim & Fabbiano 2003, 2004; Kim et al. 2006). Our procedure for correcting 



1 This V-magnitude value together with those quoted in Sect. 14.21 have been computed by Kundu & 
Whitmore (2001) assuming the Burstein & Heiles (1982) reddening. 



3.2. The LMXB X-ray Luminosity Function 




(1) 
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XLFs derived from Chandra data was developed and first applied to NGC1316 by Kim & 
Fabbiano (2003; see the Appendix therein, where the method is discussed in detail), and 
was then applied to a sample of early- type galaxies by Kim & Fabbiano (2004). A similar 
procedure had also been used to correct the XLF of the Antennae galaxies (Zezas & Fabbiano 
2002). In brief, we simulated ~ 20, 000 point sources, added them one by one to the observed 
image, and then run wavdetect to determine whether the added source was detected. With 
this procedure we correct simultaneously for incompleteness near the detection threshold, 
Eddington bias (Eddington 1913) and source confusion. In the simulations, we assumed a 
typical XLF of differential slope of f3 = 2; however, the adopted slope does not affect the 
results significantly, and neither does the assumed radial distribution of point sources (which 
we assume to follow the optical light, as observed in many elliptical galaxies, see Fabbiano 
2006). We excluded central sources with R < 10" in Tab. [31 because of the large photometric 
error caused by confusion with other overlapping sources, and by the possible presence of 
some diffuse emission; these conditions make the incompleteness corrections uncertain. With 
the above procedure, we established the 90% completeness limit (i.e., the luminosity at which 
10% of sources would not be detected inside the D25 ellipse, excluding the central 10") to 
be L x = 3 x 10 37 erg s" 1 (or F x = 4.3 x 10~ 16 erg cm" 2 s" 1 ) in the 0.3-8 keV band. 
We can reliably correct the XLF to a luminosity of ~ 2/3 of the 90% limit. The resulting 
cumulative XLF [i.e., the number of sources N(> Lx)} is shown in Fig. HI From the log N - 
log S relation derived for the ChaMP plus CDF data (Kim et al. 2007), we then estimated 
that the number of expected cosmic background sources (mostly AGN) falling within the 
D25 ellipse and brighter than the 90% completeness limit is 3, corresponding to a source 
contamination of 7%. 

To establish the functional form of the XLF, we fitted the bias-corrected differential 
XLF [eq. (1)] with single power-laws, using both the Cash and x 2 statistics in the CIAO 
Sherpa modeling and fitting application. We found that a single power-law represents well 
the XLF, with a best fit slope /3 = 2.0 ± 0.3, compatible with other power-law fits of LMXB 
XLFs (Kim & Fabbiano 2004; Kim et al. 2006) and a best fit amplitude k = SAtH- Even 
without the correction for incompleteness, that is significant only at the point corresponding 
to the faintest fluxes (Lx = 2 x 10 37 erg s _1 ), a single power law is a good fit all the way down 
to Lx = 3 x 10 37 erg s -1 . Figure H] may suggest a low luminosity break around Lx ~ 5 x 10 37 
erg s _1 . This break was suggested to be a universal feature of LMXB populations (Gilfanov 
2004), but is not observed in all cases (Kim et al. 2006). For NGC821 there is no compelling 
evidence for such a break, given the size of the error bars, and the good fit with a single 
power-law model. The combined LMXB XLF (Kim & Fabbiano 2004; Gilfanov 2004) shows 
a high luminosity break at Lx > (5.0 ± 1.6) x 10 38 erg s _1 , that may originate from the 
presence of neutron star and black hole binary populations. We cannot constrain here the 
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high luminosity shape of the XLF because of the small number of luminous sources (only 2 
sources have L x > 5 x 10 38 erg s _1 ). 

In contrast to the uniform slope, the amplitude of the XLF varies widely from galaxy 
to galaxy, reflecting the varying content of LMXBs per unit galactic luminosity. For their 
sample of early type galaxies, from integration of the XLF above Lx,min = 10 37 erg s _1 , Kim 
& Fabbiano (2004) estimate the total Lx due to LMXBs and the average ratio Lx/Lb = 
(0.9 ± 0.5) x 10 30 erg s" 1 L B l & . For NGC821 this ratio is 0.4 x 10 30 erg s" 1 L~ B l & (see 
Sect. 15.51) . lower than the mean value by la. This ratio is known to correlate with the total 
specific frequency of globular clusters, but this quantity unfortunately is not available for 
NGC821. Based on the L X /L B ratio, we would expect NGC821 to have a relatively poor 
GC population. 

4. The central region of NGC821 — resolved sources 

Figure [5] shows the central region of the image obtained from the merged data, in three 
spectral bands (0.3-1 keV, 1-2 keV and 2-4 keV) that cover the energy range where most of 
the counts are detected. From these three images in different spectral bands we also obtained 
an adaptively smoothed (using the CIAO task 'csmooth') color image (Fig. [6]), where the 
data were smoothed using scales ranging from 1 to 20 pixels (0"5 — 10"). Given the small 
extent of the region considered, no exposure correction was needed. Figure [6] shows diffuse 
emission at the center, that follows the optical shape of the galaxy: it is elongated in the 
NW-SE direction, roughly aligned with the optical major axis (see also Tab. CD, Fig. [1] and 
Fig. [H] below). Most of the diffuse emission is concentrated well within one R e (that is 43"9, 
while the size of Fig. is 55"). Its color indicates emission mostly in the soft (red) and 
medium (green) energy bands. Harder emission from a few central point-like sources is also 
visible. 

Below we concentrate on the X-ray analysis of the sources detected in the central region; 
in Sect. Owe report the analysis of the diffuse emission, for the central region and the whole 
galaxy. 

4.1. Detection and spatial properties 

We can now establish with higher significance than in F04 the morphology of the central 
X-ray emission. In particular, F04 had identified three possibly connected extended central 
emission regions (that they called SI through S3) and a point-like source to the north- 
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east of this grouping (labelled NE source), wavdetect detects all these emission regions in 
the merged data, plus another source that we call S4 (Fig. [7]). F04 had found the spatial 
distribution of photons of the S1-S3 sources to be more extended than expected for point- 
like sources seen through the Chandra mirrors, and suggested that they may be part of an 
S-shaped jet or filament, possibly connected with low-level or past nuclear activity. With the 
longer exposure, it is apparent (Fig. [7j) that the emission of S3 may be point-like, suggesting 
that this source is a luminous LMXB, but the other sources are more extended. 

In order to establish the spatial properties of S1-S4, following F04 we compared the spa- 
tial distribution of their counts with that of the on-axis image of the quasar GB 1508+5714, 
a hard point-like source (Siemiginowska et al. 2003). GB 1508+5714 gives a good repre- 
sentation of the Chandra ACIS-S PSF for our analysis, since it lies at the same distance 
from the aimpoint as S1-S4 and is similarly hard (see Sect. 14.31 below). The image of 
GB 1508+5714 contains 5,300 counts within 2" of the centroid of the count distribution; 
the ratio of counts within the 1"— 2" annulus to those in the central circle of 1" radius is 
Ratio{PSF)=0.043±0.001 (la). The isolated NE source in this central field (Fig. has a 
total of 173 source counts and Ratio(NE)=0. 123+0.083, a value consistent with that of our 
reference quasar, within the la error. 

The total number of source counts for S1-S4 are 184, 250, 326 and 84 respectively; the 
analogous ratios for the background-subtracted counts, around the centroids determined by 
wavdetect, are Ratio{Sl)= 0.59+0.11, Ratio{S2)=0. 77+0.10, Ratio{S 3) =0.20±0. 06, and 
i?atzo(S'4)=0.80±0.15. Therefore the spatial distribution of the counts from S1-S4 is ex- 
tended, when we compare it with that of the quasar. If, more conservatively, we use as 
a benchmark the NE source, we find that S3 is point-like, while SI, S2 and S4 are all 
significantly more extended than the NE source (at 3.5, 4.6 and 4a respectively). Al- 
though the significance of the NE source is such that its Ratio is less well determined 
than that of the quasar, this source would be affected by similar small amounts of smear- 
ing in the merging process as the S1-S4 sources. [Note that in F04 the NE source had a 
Ratio(NE)=0. 057+0. 054, consistent within the errors with the value we find here, although 
the nominal value is smaller; this justifies using NE as a benchmark.] 

In conclusion, excluding source S3, the counts of SI, S2 and S4 cannot come entirely 
from a point source. While the previous Chandra observation was suggestive of a central 
emission feature intrinsically elongated, possibly mostly due to diffuse emission (F04), the 
present imaging analysis suggests a cluster of extended sources. With these data it remains 
possible that each source is truly extended or produced by several point-like components, or 
alternatively due to a point-like component embedded in truly diffuse emission. 
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4.2. Position of the nucleus - S2 

From the Two Micron All Sky Survey (2MASS, Skrutskie et al. 2006), the center 
of NGC821 is at RA=02 A 08 m 21.14, Dec=+10° 59' 4T/7 (J2000), with an uncertainty of 
l'/25 (at 95% confidence), as reported in the NASA/IPAC Extragalactic Database (NED). 
Considering also the Chandra aspect nncertaintjl (the 90% uncertainty circle has a radius 
of < 0'.'6, and the 68% circle of < 0'.'25), this position agrees with that of the extended source 
S2 determined by wavdetect, that is RA=02 h 08 m 21110, Dec=+10° 59' 41'.'6, with an offset 
of just 0'.'6. 

In addition, we compared the Chandra and HST positions by using the archival WFPC2 
F555W and F814W filter images (also used by Kundu & Whitmore 2001, Sect. 13.11) . First we 
improved the astrometry of the WFPC2 images, re-fitting the coordinate grid to all sources 
with positional errors < 0.'3 in the 2MASS All-Sky Catalog of Point Sources (Cutri et al. 
2003) and the USNO-B1.0 Catalog (Monet et al. 2003). We found 5 coincidences between 
the WFPC2 and the USNO-B1.0 sources, two of which are also 2MASS sources; with these, 
we improved the HST absolute astrometry to < 0'.' 2 in the nuclear region. We then compared 
the recalibrated WFPC2 images with the Chandra image, and identified eight sources with 
optical/X-ray coincidences, not including the nuclear source S2 (Fig. [8]). Of these, seven lie 
inside the D25 ellipse and have been already picked up in Sect. l3.ll as associated with globular 
clusters, with an offset smaller than the uncertainty in the X-ray position (Tab. [3D; however, 
for the present astrometric analysis, we excluded one of them (CXOU J020817. 7+105907) 
located almost at a chip edge. The remaining six coincidences are shown with yellow circles 
in Figure [8j They have nearly pointlike appearance, as expected for GCs at the distance of 
NGC821, and their V-magnitudes range from 20.9 to 22.4, i.e., they are among the brightest 
GCs of NGC821 (Kundu & Whitmore 2001, as also discussed in Sect. 13. ip . The eighth source 
with optical/X-ray coincidence is much brighter (by almost one magnitude in the V-band) 
than the brightest GC, and lies outside the D25 ellipse (it is located at RA=02 /l 08 m 17172, 
Dec=+ll° 00' 3i'.'22); it is unresolved and could be a background AGN, as there is a higher 
chance of a background source outside the optical image. 

There is no systematic shift or rotation between the revised HST and Chandra posi- 
tions. For each of the seven coincidences, the displacement between the optical and X-ray 
positions is < / .'22, and the root-mean-square displacement for the whole sample is a = 0'.'14. 
The relative uncertainty between HST and Chandra is therefore good to 0'.'14, which sug- 
gests that the Chandra pointing was in fact very accurate. In the recalibrated WFPC2 image 
the optical nucleus is located at RA=02 /l 08 m 21 .13, Dec=+10° 59' 41'.'8, that is within 0'.'2 



2 The Chandra absolute astrometric accuracy is discussed at 



http://cxc.harvard.edu/cal/ASPECT/celmon/ 
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of the 2MASS position. As shown by Fig. [9], the WFPC2 nucleus is (X'5 off the wavdetect 
center of S2, a displacement larger than the average ft' 14 offset, which can be attributed to 
the fact that S2 is an extended source. In conclusion, within the accuracy of the best relative 
astrometry available, S2 is coincident with the optical center. 

4.3. Spectral analysis of sources SI— S4 

We extracted the spectra of sources S1-S4 from the merged event file, using the extrac- 
tion regions shown in Fig. [7J the background spectrum was derived from a source-free region 
outside the optical body of the galaxy but still inside the ACIS-S3 chip. We also generated in- 
strumental response files, weighted by the exposure time of each individual observation. Prior 
to fitting, we binned the spectra in order to have at least 20 total counts per energy bin, with 
the exception of the weakest source (S4) that was rebinned at 10 counts per energy biro. We 
then fitted the spectra with an absorbed power-law model [XSPEC model: wabs(powerlaw)} 
and an optically thin thermal model corrected for absorption [XSPEC model: wabs(apec), 
with abundance Z=0.5 in solar units, and the solar abundance values of Anders & Grevesse 
1989]; the absorption was let to be a free parameter. The low number of counts did not 
allow for fits with composite spectral models. The results are summarized in Table [51 

The spectral analysis indicates a hard spectral shape (thermal emission of kT > 3 keV 
or power law emission with T = 1.4 — 2.3), with little or no intrinsic absorption. The resulting 
column densities Nh are consistent with the Galactic value (6.2 x 10 20 cm -2 ) for most of the 
fits in Table [5j if Njj is fixed to be Galactic, the best fit spectral parameters are substantially 
unchanged within the uncertainties. The spectra of S1-S4 along with their best fit power 
law models are shown in Fig. [TUJ The spectral properties of S1-S4 are overall consistent 
with those typical of LMXBs (see, e.g., Sect. I3.ip . even though these sources look extended 
(Sect. I4.ip . The source at the nucleus (S2) seems to be harder than the other sources. 

Table lists the source luminosities calculated from the best fit parameters. S2 has a 
luminosity of about 6xl0 38 erg s -1 ; since it is extended, a nuclear point source, if contributes 
to the emission of S2, will have a smaller luminosity. Since we have now established that S2 
is coincident with the nucleus (Sect. 14.2)) . it is reasonable that a fraction of the S2 emission 
originates from the MBH. However, an estimate of this fraction would not be unique, because 
it will require some model assumptions (e.g., one could model S2 as a point source plus some 
extended emission in various ways). F04 derived a 3 a upper limit for the nuclear luminosity 



3 When grouping with different criteria, or using the Cash statistics instead of the % 2 one, the spectral 
results for this source do not change significantly. 



13 



of Lx < 4.6 x 10 38 erg s" 1 in the 0.3-10 keV band; they assumed that all of the detected 
emission in the central S-shaped feature (Sect. 14. ip was truly extended, and calculated the 
limit from the largest emission found in a 1" x 1" sliding cell moved across this feature. 
In order to avoid model dependent results, we adopt as upper limit to the MBH emission 
the value derived with the same procedure used in F04, with the 1" x 1" cell slid over and 
around S2. This gives a 3cr upper limit of Lx < 2.8 x 10 38 erg s _1 in the 0.3-8 keV band, 
and L x < 1.8 x 10 38 erg s" 1 in the 2-10 keV band. 

5. The diffuse emission and the hot ISM contribution 

Estimating the amount of gaseous emission in NGC821 is important, because this hot 
gas could be a source of fuel for the nuclear MBH (see Sect. [Q). The ACIS-S adaptively 
smoothed mapped color image of the NGC 821 field (Fig. Ej) shows diffuse emission at the 
center, roughly aligned with the major axis of the galaxy. However, this diffuse emission is 
not by itself proof of the presence of a hot ISM, since an important contaminant in a galaxy 
like NGC821, that has a low hot gas content (Sect. [TJ, is the undetected portion of the 
LMXB population (see for example the case of NGC1316, Kim & Fabbiano 2003); another 
contaminant could be the emission of normal stars (first discussed for this type of galaxies 
in Pellegrini & Fabbiano 1994). 

In this Section we report the results of the analysis of this diffuse emission. First, we 
analyze the merged data to estimate its spatial and spectral properties (Sections 15.11 and 
15. 2p . Since our deep data allow for an accurate characterization of the LMXB population 
(Section EJ), we then use these results to estimate how much of the diffuse emission could 
come from unresolved LMXBs below the detection threshold, comparing: (1) the shape of 
the radial profiles observed for the diffuse emission, the resolved point sources and the stellar 
optical light (Sect. 15.31) ; (2) the radial profile of the soft emission contributed by LMXBs 
and the total observed soft profile (Sect. I3T4T) ; (3) the total luminosity of the diffuse emission 
with that expected from unresolved LMXBs, using the XLF derived in Sect. 13.21 (Sect. 15. 5p . 
Finally, in Sect. 15.61 we derive the expected contribution from other types of stellar sources. 

5.1. Radial profile 

After having excluded all the detected sources, we derived a radial profile of all the 
remaining counts, in the 0.3-6 keV energy band; at energies > 6 keV the background con- 
tribution becomes dominant and the addition of these photons would only increase the 
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errorbars. The profile was centered on the nucleus, and included 30 circular annuli (each 
~ 5" wide). The resulting radial profile is shown in Figure [TTJ where the presence of emission 
is evident as far as ~ 20" — 30" from the nucleus (~ 2.3 — 3.5 kpc at the galaxy distance). 
The flattening at larger radii is consistent with the expected background level. Therefore, 
diffuse emission is detected out to a radius that lies well within the optical extent of the 
galaxy, even within one effective radius (Tab. [T|). 

5.2. Spectral analysis of the diffuse emission 

We analyzed the spectra of three circles centered on the galactic center, of progressively 
larger radius: 5" and 10", to investigate the gas presence in the circumnuclear region, and 30", 
that is the outermost radius at which diffuse emission is detected. Prior to their extraction, 
we removed regions including detected sources. The background was estimated from a source- 
free circular area of 50" radius, located outside the galaxy. We fitted the spectra with single 
component models, such as a power law (XSPEC model wabs * pow) as expected from the 
integrated LMXB emission, and a thermal APEC model (XSPEC model wabs*apec), to 
account for the emission of a hot gas component; we then used a composite power-law and 
APEC model (wabs(apec+pow)) . The results are summarized in Tab. El errors quoted below 
give the 68% confidence interval for one interesting parameter. 

The spectrum of a central circle of 5" radius is well fitted by a single absorbed power 
law model, with a photon index T = 1.77 ± 0.17 and Nh consistent with the Galactic value, 
which is in good agreement with the spectral parameters of LMXBs (Sect. |3~TI) . The thermal 
model also gives an acceptable fit, provided that the temperature is high (kT ~ 5 keV at the 
best fit), which again is what is expected for LMXBs (Fabbiano 2006). Next we attempted 
the wabs(apec+pow) spectral model, to probe for the presence of a hot gas component, in 
addition to the LMXB emission that is clearly dominant. Nh was fixed at the Galactic value 
and the abundance to the solar value, that is an average for the stellar population in the 
central galactic region (Proctor et al. 2005). At the best fit, T = 1.61+^ and kT = 0.21+°;^ 
keV; however, this fit does not represent a statistically significant improvement with respect 
to the simple power law model, as was established via a "calibration" of the F-test through 
simulations, following the prescriptions^ discussed in Protassov et al. (2002; their Sect. 5.2). 



4 The basic steps of the procedure consist of fitting the observed data with the simple (pow) and composite 
(apec+pow) models, getting the F-statistic for these fits, simulating a large number spectra for the simple best 
fit pow model with XSPEC, fitting each of the simulated data sets with both the simple and the composite 
model and recording the F-statistic, getting the distribution of the simulated F-statistic to compute a p- value 
(Protassov et al. 2002) that establishes whether the addition of the thermal spectral component is significant 
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The analysis of the spectrum of a central circle of 10" radius gave similar results (Tab. [6]). 
This spectrum is well described by a power law consistent with the average spectrum of 
LMXBs (r = 1.81±0.13 and Galactic Nh), and again a thermal component is not statistically 
required. 

We then analyzed the spectrum of the diffuse emission from the whole galaxy (i.e., from 
within a radius of 30"). Again, a simple power law model gives an acceptable fit, with 
T = 1.74 ± 0.14. If a thermal component is added, with solar abundance and Galactic Nh, 
then T = 1.59^Q °y and kT = 0.59±0.18 keV; however, this thermal component is not required 
statistically. The 68% confidence upper limit on its luminosity is L(0.3-8 keV)< 1.34 x 10 38 
erg s -1 . Based on the shallower pointings of 2002, David et al. (2006) estimated a 90% 
confidence upper limit on L(0.5-2 keV) of 7.0 x 10 38 erg s _1 , for a thermal component. 

In the analysis above the presence of a small amount of thermal emission confined to 
the central region might be overwhelmed by the emission of unresolved LMXBs along the 
line of sight. Therefore we also attempted the deprojection of the spectral data by using the 
technique implemented within XSPEC, where the spectra extracted from annuli centered 
on the nucleus are compared with the spectra expected from the superposition along the 
line of sight of the emission coming from the corresponding spherical shells [as successfully 
done for the Sombrero bulge, Pellegrini et al. (2003b)]. The central region was therefore 
divided in an inner circle of 5" radius, a surrounding annulus with outer radius of 10", and 
another surrounding annulus of outer radius of 20". The spectra of each of these regions were 
jointly fitted with the spectral models used above. Unfortunately, the outermost annulus 
did not have enough counts to add meaningful spectral information, therefore only the two 
inner regions were deprojected. The simple power law model (projct*wabs*powerlaw) , with 
the photon index constrained to be the same in the two regions, gives an acceptable fit 
for T = 1.68i{u5 and Nh = 2.511 jj x !0 20 cm" 2 (inner region) and < 4.3 x 10 20 cm" 2 
(outer region). The addition of a thermal component [XSPEC model projct*wabs(apec+pow), 
with solar abundance and Galactic Nh] gives a best fit T = 1.73l£^, and temperatures of 
kT = 0.17lo'o7 keV for the outer annulus, and kT = 0.031qo3 keV for the inner circle. 
The latter value is too small to be meaningful, therefore the fit was repeated with a fixed 
T = 1.7 (typical of LMXBs, Sect. I3.ip . The outer temperature remained unchanged and a 
more meaningful temperature was obtained for the inner circle (kT = 0.18 keV), but it was 
basically unconstrained. Again, the addition of the thermal component was not statistically 
significant. At the best fit, the luminosity of the thermal component is ~ 5% and ~ 13% 
that of the power law, respectively for the inner and outer regions. 



or not. 
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In summary, the analysis of the projected regions and of the central deprojected regions 
never requires the presence of a soft thermal component at a statistically significant level. 
The luminosity of a soft thermal component inserted in the fit is much lower than that of 
the power law (< 1/10 at the best fit) for the spectra of the 5", 10" and 30" circles. 

5.3. Radial profiles of detected sources, diffuse emission and galactic optical 

emission 

Figure [T2] shows the comparison of the radial profiles of three quantities: (1) the 
background-subtracted diffuse emission in the 0.3-6 keV band, (2) the number density of 
resolved point sources, (3) the R-band emission of NGC821. The latter, that gives the dis- 
tribution of the stellar surface brightness, is well described by a de Vaucouleurs (1948) law 
I{R) oc exp {-7.669 * [(R/43^9) 1 ^ 4 - 1]}, where R is in arcseconds (Soria et al. 2006b). 

All three profiles follow the same radial trend. The agreement between the radial profile 
of resolved point sources and that of the stellar light is not surprising (see, e.g., Kim & 
Fabbiano 2003, Fabbiano 2006). The similarity of these profiles with that of the diffuse 
emission provides instead support to (or at least is consistent with) the idea that also the 
diffuse emission mostly comes from undetected LMXBs. 

5.4. Simulated (LMXB) and observed soft radial profiles 

If the diffuse emission is mostly due to unresolved LMXBs, as suggested by the analysis 
of the previous two Sections 15.21 and I5.3[ we would expect its spectrum to have a power- law 
shape typical of the LMXB population, at all radii. Any deviation is suggestive of a localized 
additional hot ISM component. 

In order to further check the lack of a soft thermal component indicated by the previous 
analysis at all radii, we made the following test. We derived from observations the radial 
profile of the diffuse hard emission in the 1.5-6 keV band; in the hypothesis that it is 
produced by a distributed source whose spectrum is a power law of T = 1.7 and Galactic 
Nh, we then used it to derive a "simulated" radial profile in the 0.3-1.5 keV band. The chosen 
spectrum represents the emission of unresolved LMXBs, in accordance with the results of 
Sects. 13.11 In Figure [TBI we compare the simulated soft profile with the observed profile of the 
diffuse emission in the same energy band. The agreement between the two profiles is good, 
consistent with the idea that the observed soft emission is due substantially to unresolved 
binaries. Only within the central circle of 10" radius an additional contribution from a soft 
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source could be present, but is not revealed by the spectral analysis (Sect. 15 .2j) . 

5.5. The unresolved emission implied by the XLF 

The emission from undetected LMXBs can be recovered from the XLF derived in 
Sect. 13.21 specifically for NGC821, as the difference between the total luminosity expected 
from LMXBs based on the XLF, Lx (total), and the luminosity due to all detected sources, 
Lx (detected). In the 0.3-8 keV band, Lx (detected) = 6.5 x 10 39 erg s _1 for the sources falling 
within the D25 ellipse. For the same region, Lx (total) is not just given by an integration of 
the XLF, since the XLF was derived for the D25 ellipse with a central circle of 10" radius 
excluded (Sect. 13. 2j) . Lx (total) for the whole D25 ellipse is then obtained rescaling the lu- 
minosity resulting^ from the XLF by the ratio between the luminosity of detected sources in 
the two regions (i.e., the D25 ellipse and this ellipse without the central circle). This proce- 
dure relies on the assumptions that (1) the XLF for R < 10" and R > 10" is the same, and 
(2) Lx (total) / L x (detected) is also the same in these two regions. While (1) is reasonable, 
(2) is less valid, since more sources would be hidden in the inner region. Therefore, both 
Lx (total) and Lx (undetected) =Lx (total) - Lx (detected) are actually lower limits. For the 
D25 ellipse, Lx (total) is 7.6 x 10 39 erg s" 1 , and L x (undetected) = 1.1 x 10 39 erg s" 1 . 

It is interesting now to compare the unresolved LMXBs emission with the luminosity of 
the diffuse emission derived from the spectrum of a central circle of 30" radius (i.e., 1.4xl0 39 
erg s -1 , Sect. 15.21) . Therefore we repeated the procedure above to calculate Lx (undetected) 
for this circle, and it turned out that Lx (undetected) =7.4 x 10 38 erg s _1 . This value is slightly 
lower than the spectral luminosity of the diffuse emission; however, as noted above, it is likely 
an underestimate of the true value. Also, considering the uncertainty in the normalization 
of the XLF (k = Sect. 153]) . L x (undetected) can be as large as 2.7 x 10 39 erg s~\ 

In summary, for the same circle of R = 30", the spectral luminosity of the diffuse 
emission and the XLF-based unresolved luminosity are close, and they agree within the 
uncertainties. 



5 This integration assumes that the differential XLF has the same slope (~ 2) down to Lx = fO 37 erg s _1 , 
and is done over the range Lx — (10 37 — 10 39 ) erg s^ 1 . Given the slope value, sources of luminosities below 
10 37 erg s _1 do not contribute significantly to the result; for example, decreasing the lower boundary by a 
factor of two would increase the resulting total luminosity by ~ 15%. In addition, a low luminosity break in 
the XLF is also expected (e.g., Kim et al. 2006). 
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5.6. Other stellar sources 

Since we are constraining the origin of the diffuse emission in NGC821 better than it 
has been possible so far in previous analyses of low Lx/ Lb ellipticals, we check also for the 
amount of diffuse emission that can be accounted for by stellar sources other than LMXBs. 
In the old stellar population of an elliptical galaxy these include coronae of late type main 
sequence stars, RS Canum Venaticorum (RS CVn) systems and supersoft X-ray sources 
(Pellegrini & Fabbiano 1994). The latter are characterized by black body emission with 
effective temperatures of 15-80 eV, X-ray luminosities up to fewxlO 38 erg s _1 and are more 
frequent in late type and irregular galaxies (Kahabka & van den Heuvel 1997; Di Stefano & 
Kong 2004). NGC4697 is the only elliptical where they have been identified so far (Sarazin 
et al. 2001). In the X-ray color-color plot of NGC821 (Fig. [2]), supersoft sources would fall 
in the upper right corner, and there is no such source with > 30 counts; one or two sources 
with less than 30 counts have a high C32 color, but with a large error. The contribution of 
supersoft sources will not be considered further here. 

Stellar coronae of main sequence stars are all sources of thermal X-ray emission {Lx ~ 
10 26 — 10 30 erg s _1 ) indicative of coronal plasmas at temperatures of ~ 10 6 — 10 7 K (e.g., 
Schmitt et al. 1990). These luminosities are faint in comparison to accretion powered 
stellar sources, but main sequence stars of spectral types G, K and M are present in a very 
large number (hereafter respectively Nq, Nr and Nm)- These N{ can be recovered from 
the initial stellar mass function (IMF) ip(M) = AM' 2 - 35 (Salpeter 1955), where the scale- 
factor A = 1.67Lb, with Lb in Lg jQ , for a 12 Gyrs old stellar population with a 0.5 solar 
metallicity (Maraston 2005), as suitable for NGC821 (Proctor et al. 2005). Then N t = 
A j^ su ^ M~ 2 - 35 dM, where M in f and M sup are respectively O.7M and the main sequence 
turn-off mass of O.9M (as appropriate for NGC821, Maraston 2005) for G stars; 0.5 and 
O.7M for K stars; and 0.1 and 0.5M Q for M stars. Integration gives Nq = 0.6L#, Nk = 
I.ILb and Nm = 24.5Lb, therefore the collective M dwarf X-ray emission is by far the most 
important, because the average X-ray luminosity of a G, K and M star (hereafter Lx,m) are 
comparable (Kuntz & Snowden 2001). Taking Lx,m J$ 1-4 x 10 27 erg s _1 in the 0.3-8 keV 
bancu, the collective M dwarfs emission in NGC821 is then Nm x Lx,m ^ 6.4 x 10 38 erg s _1 . 
From a central circle of 30" radius, for the de Vaucouleurs optical profile of Sect. I5.3[ the M 
dwarfs emission is < 2.5 x 10 38 erg s" 1 , that is < 18% of the luminosity derived from the 
spectrum of the diffuse emission from the same region (Sect. 15. 2p . However, recent claims 



6 The < sign comes from assuming thermal emission of kT = 1 keV when converting Lx,m in the ROSAT 
band (Kuntz & Snowden 2001) to the 0.3-8 keV band. The proper spectral description consists of two thermal 
components: one with T ~ 2 — 4 x 10 6 K and the other with T ~ 10 7 K; the latter usually has larger (and 
often much larger) emission measure (Giampapa et al. 1996). 
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favor an IMF natter than the Salpeter one at the low mass end, as for the Kroupa (2001) 
IMF, which is oc M~ L3 for 0.08 < M/M < 0.5. In this case the scale-factor A = l.69L B 
(Maraston 2005), the integration gives Nm = 4.3L^, and the collective M dwarfs emission 
becomes < 1.1 x 10 38 erg s _1 for the whole galaxy, and < 0.5 x 10 38 erg s _1 for a circle of 
30" radius. Therefore, it can account for < 3% of the luminosity of the diffuse emission. 

RS CVn systems are chromospherically active objects, consisting of a G or K giant or 
subgiant, with a late type main sequence or subgiant companion (Linsky 1984). They are the 
most X-ray luminous late type stars (Lx ~ 10 29 — 10 3L5 erg s _1 ), and their spectra can be 
well modeled by a thermal plasma with two temperatures (of average kT =1.3 keV and 0.18 
keV, Dempsey et al. 1997). From Dempsey et al.'s Fig. 3, the average 0.3-8 keV luminosity 
of a system id_l < 3 x 10 30 erg s -1 . The number of RS CVn systems is a fraction of the number 
of giants and subgiants (respectively N giant and N subg ) expected to be present in NGC821, 
which is^| N su b g + N giant ~ 0.04L#. Taking half this number for those giants or subgiants 
that are also in binary systems, and another factor (~ 0.2) for those that become RS CVn, 
then the collective contribution is < 0.1(N subg + N giant ) x 3 x 10 30 erg s _1 = 1.2 x 10 28 L S 
erg s . This is < 2.2 x 10 38 erg s _1 for the whole galaxy, and < 8.9 x 10 37 erg s _1 from 
a circle of 30" radius, that is < 6% of the total diffuse emission within the same region. A 
somewhat higher percentage (~ 15%) is derived in the 2-10 keV band when adopting the 
emissivity per unit stellar mass of these sources in the Galactic plane (Sazonov et al. 2006), 
and assuming an average stellar mass-to-light ratio in the B-band of ~ 6.5 (Sect. 16. II) . This 
higher percentage could be in part the result of a higher emissivity in the Milky Way, due 
to younger stellar ages than estimated for NGC821 (see above). 

In conclusion, the collective emission of stellar sources other than LMXBs from the whole 
galaxy is < 3.3 x 10 38 erg s -1 ; as expected, this is much smaller than the collective LMXB's 
emission (7.6 xlO 39 erg s" 1 , Sect. but it is even smaller than its unresolved fraction 

[i.e., Lx (undetected) = 1.1 x 10 39 erg s _1 , Sect. 15.5] . Also, these sources can contribute up to 
~ 10% of the diffuse emission detected within a radius of 30" (Sect. 15.21) . 



7 Again the < sign comes from assuming that most of the emission is due to the higher temperature 
component. 

& Ngi an t and N su b g are the product of the time spent by stars in such evolutionary phases by the specific 
evolutionary flux for a 12 Gyr old stellar population (which is fairly independent of the IMF; Maraston 2005) 
and by Lb- Low mass metal rich stars spend 7 x 10 s yrs in the G/K-subgiant phase, and 5 x 10 8 Gyr on 
the red giant branch (e.g., Renzini 1989). 
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5.7. Conclusions on hot gas presence 

All the different lines of investigation undertaken in this Section consistently indicate a 
largely dominating LMXB contribution to the diffuse X-ray emission. Its spectrum does not 
require soft thermal emission, that is limited to < 1.34 x 10 38 erg s^ 1 (< 10% of the total 
diffuse luminosity) at 68% confidence. The observed radial profiles of the stellar light, the 
resolved point sources and the diffuse emission agree within the errors, suggesting that the 
diffuse emission has a stellar origin. The observed radial profile in the soft band is consistent 
with that expected from the sources (the LMXB population) producing the hard emission, 
except perhaps for the innermost (R < 10") region. Finally, the luminosity of undetected 
sources derived from the bias-corrected XLF agrees with the spectral luminosity of the diffuse 
emission, within the uncertainties. A fraction of the diffuse emission (< 10%) can also come 
from stellar sources other than LMXBs. 



6. Discussion 

In the previous Sections we have reported the results of the analysis of a deep Chandra 
pointing (for a total exposure of nearly 230 ks) aimed at detecting the nuclear emission 
and setting stringent constraints on the presence of hot gas both on the circumnuclear and 
the galactic scale. The deep Chandra image revealed a significant portion of the LMXB 
population, for which the X-ray Luminosity Function was derived. This, together with a 
spectral and imaging analysis of the diffuse emission, severely constrained the presence of a 
hot ISM (Sect. [3]). At the galactic center, in addition to diffuse emission (which we can explain 
as predominantly due to unresolved LMXBs), we find four sources of Lx = (2 — 9) x 10 38 
erg s _1 . We have been able to establish that one of these sources (S2) is coincident with the 
galactic center. Three of these source, including S2, are extended, with one of them as long 
as 5" (SI; Fig. [7]). These sources could be truly extended, or made of one or more point- 
like components, possibly also embedded in truly diffuse emission. Their spectral shape is 
consistent with that of LMXBs, but in NGC821 there are only 2 other sources as bright as 
these, outside the central circle of 10" radius (Fig. H|). This supports the idea that more 
than one LMXB is contributing to the emission of each source. An alternative possibility 
is that they belong to a continuous complex feature, made of LMXBs (that may also lie 
there due to a projection effect), a nucleus (possibly with associated extended emission) 
and an outflow/jet. SI, in particular, has a linear structure suggestive of a jet (Fig. [7]). 
Our VLA observations (Sect. [I]) help better address these possibilities, and are presented in 
the companion paper (Pellegrini et al. 2007), together with a discussion on whether/how 
accretion proceeds in this nucleus, based on all the available observational evidence. We 



-21 - 



estimated here a 3cr upper limit to the 0.3-8 keV emission of a point-like source associated 
with the MBH of < 2.8 x 10 38 erg s _1 , which makes this one of the quietest MBH studied 
with Chandra (e.g., Pellegrini 2005b), with L x /L Edd < 2.5 x 10~ 8 . 

In the following we discuss the possibility that there may be undetected hot gas in 
NGC 821 (Sect. 16.11) . Hidden hot gas on the nuclear scale could be a source of fuel for the 
MBH; furthermore, we check whether the lack of detection on the galactic scale agrees with 
expectations for the evolution of the ISM in this elliptical. 

6.1. Is there hot undetected gas? 

The present analysis does not detect hot gas available for accretion at the nucleus; 
coupled with the absence of interstellar medium observed at other wavelenghts (Sect. [Tj), the 
lack of fuel for the MBH could be the simplest explanation for the very low levels to which the 
nuclear emission has been constrained. However, an aging stellar population continuously 
returns gas to the ISM, via stellar mass losses (e.g., Ciotti et al. 1991, David et al. 1991); 
as a mimimum, the circumnuclear region should be replenished with this fuel. Below we 
discuss the possible presence of this ISM. 

The lack of detection of hot gas on the galactic scale can be expected from simple 
energetic considerations, comparing the energy made available per unit time by type la 
supernova explosions (Lsn) with the energy required to steadily extract from the galactic 
potential well the stellar mass losses produced per unit time (L grav ; Ciotti et al. 1991), 
for the whole galaxy lifetime. Lsn comes directly from the observed SNIa's explosion rate 
in ellipticals (Cappellaro et al. 1999) rescaled for Lb in Tab. (TJ this rate was higher in 
the past (e.g., Greggio 2005). L grav derives from the stellar mass loss rate and the galactic 
mass profile, that can be well modelled by the superposition of two Hernquist (1990) density 
profiled, one for the stars with total mass M*, and one for the dark matter with mass M^. 
This model is then tailored onto NGC821 taking the observed L#, R e , central stellar velocity 
dispersion (Tab. [T]) and M^/M* ratio. Anisotropic Jeans models reproducing the observed 
velocity dispersion profiles of stars and planetary nebulae out to 5R e constrain the range 
of (Mh + M*)/Lb to be 13-17 (Romanowsky et al. 2003). For a plausible stellar mass-to- 



9 This model gives a very good approximation of the de Vaucouleurs (1948) law, that fits the light profile 
of NGC821 (Sect. 15.3)) . The radial distribution of the dark haloes of ellipticals is not well constrained by 
the observations. Theoretical arguments favor a peaked profile (Ciotti & Pellegrini 1992; Evans & Collett 
1997). High resolution numerical simulations (Navarro, Frenk, & White 1996) produce a central density 
distribution equal to that of the Hernquist model. 
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light ratio M*/Lg (Gerhard et al. 2001, Napolitano et al. 2005) one derives M h /M* ~ 1.2. 
Finally, the stellar mass loss rate is distributed as the stars, and derives by summing the 
mass ejected by each star as a function of its mass (Renzini & Ciotti 1993) for a Kroupa 
IMF, scaled for NGC821 as described in Sect. I5.6( this produces a total mass loss rate of 
0.25 M Q yr" 1 after 12 Gyrs, in good agreement with the values estimated for ellipticals from 
ISO observations (Athey et al. 2002). The resulting Lsn ~ (3 — 4)L grav during the galaxy 
lifetime, with L SN ~ 4L grav at present; therefore type-la SNe have always provided the 
heating to drive the stellar mass losses in a galactic wind (neglecting the effect of radiative 
cooling that is not expected to be important for very low gas densities typical of a wind). 

However, this is a global energetic calculation, and the galaxy may host regions with 
L SN > L grav and regions with L SN < L grav . For centrally peaked mass distributions like 
those described above, numerical simulations of hot gas evolution showed that this is indeed 
the most frequent case (Pellegrini & Ciotti 1998): the ISM is driven out of the outer galactic 
regions in a wind, while it is inflowing within a stagnation radius that can be largely differ- 
ent from galaxy to galaxy. To explore whether this situation applies to NGC 821, we ran 
hydrodynamical simulations specific for the galaxy model described above, plus the addi- 
tional gravitational attraction contributed by its central MBH (see Pellegrini & Ciotti 2006 
for more details on the numerical code and the time-evolving input quantities). The central 
grid spacing was set to 5 pc to allow for a better sampling of the inner regions and the gas 
flow evolution was followed for 12 Gyrs, an age comparable with the stellar age of NGC 821 
(Proctor et al. 2005). The flow kept in a partial wind all the time, and at the end its total 
Lx ~fewxl0 36 erg s _1 in the 0.3-8 keV band, well below our observational limit on the hot 
gas. The bulk of the hot gas was outflowing, and directed towards the center from within a 
radius of ~ 25 pc. At the smallest radius at which the flow is well resolved (10 pc) the mass 
inflow rate is M in fewx 10 _5 M Q yr _1 ; the gas density is n e ~ 10~ 3 cm -3 , one third of the 
density value calculated from the emission measure of the "best fit" thermal component for 
the central circle, in the deprojected analysis of the diffuse emission (this component is not 
statistically required, Sect. 15.21) . Note however that the inflowing region is very small and 
fully included within the extent of the S2 source (Sect. 14.11) . By varying the M^/M* ratio, 
the SNIa's rate, the age of the galaxy and the stellar mass loss rate within the limits imposed 
by observational uncertainties, the value of M in keeps within the range (2 — 7) x lO _5 M yr _1 . 
The implications for the nuclear emission coming from this estimate of M, n are discussed in 
Pellegrini et al. (2007). Here we just note that Mj„ ranges from 10~ 3 to fewxlO~ 5 M yr _1 
in the past ~ 10 Gyrs, and its integration over this time gives a total mass accreted at the 
center of ~ 2 x 1O 6 M , that is 2.3 x 10~ 2 of the observed MBH mass (Tab. HJ. Therefore, 
accretion of stellar mass losses was not effective in building this MBH mass. 

Finally, we comment on the fact that, by using Chandra observations, the mass accretion 
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rate of galactic nuclei is customarily estimated from the analytic formula of Bondi (1952), 
valid for spherically symmetric accretion from a nonrotating polytropic gas with given density 
and temperature at infinity (e.g., Loewenstein et al. 2001, Soria et al. 2006a). Infinity is 
replaced with an accretion radius r acc = 2GMbh / c 2 s (Frank et al. 2002), where the sound 
speed c s oc \JkT is calculated as close as possible to the MBH. Even in those cases where 
kT can be observed close to the galactic center, so that one can be confident to derive an 
approximate estimate of r acc , there are additional ingredients that must be considered in the 
estimate of the mass accretion rate, which were not included in the Bondi (1952) treatment: 

1) the presence of mass and energy sources (i.e., stellar mass losses and SNIa's heating); 

2) the presence of cooling; 3) the fact that r acc is not a true "infinity" point, since the gas 
here experiences a pressure gradient. The simulations described here take into account all 
these aspects, and should give a more reliable mass accretion rate than the Bondi theory. 
However, they have some limits too. At these very low M in values, and very small inflowing 
regions, the detailed shape of the galactic mass profile becomes important. In the case of 
NGC821, the observed profile is steeper than modelled here at the very center (Gebhardt et 
al. 2003), which should produce a larger M in . The true accretion rate could be somewhat 
higher than M in also if even the stellar mass losses within the innermost radius resolved 
by the simulations are to be accreted. But further pursuing the flow behavior closer to the 
MBH with the simulations used here runs into a problem: the discrete nature of the stellar 
distribution becomes important, since the accretion time (~ 3 x 10 5 yrs from 10 pc, in the 
simulations) becomes comparable to (or lower than) the time required for the stellar mass 
losses to mix with the bulk flow (Mathews 1990), and the time elapsing between one SNIa 
event and the next. 



7. Summary 

We have observed the nearby, inactive elliptical galaxy NGC821, known to host a central 
MBH, with deep Chandra pointings, in order to put strong constraints on its nuclear emission 
and the presence of gas available for accretion. Our results can be summarized as follows: 

1. We detect 41 sources within the optical size (D25) of the galaxy; excluding 4 particular 
sources at the center, these represent the LMXB population of NGC 821. Their X-ray 
colors, and the spectral analysis of the three brightest ones, are consistent with spectral 
shapes described by power laws of photon index in the range T = 1.5 — 2.0, without 
significant intrinsic absorption. There is no clear evidence for supersoft sources. 

2. In the overlap region between D25 and the field of view of previously taken HST 
WFPC2 images, there are six clear associations of LMXBs with globular clusters, plus 
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two additional marginal matches. These LMXBs reside in the brightest GCs of the 
galaxy. Excluding the two marginal matches, ~ 30% of the X-ray sources lie in GCs, 
and ~ 6.6% of the GCs have a LMXB, broadly in agreement with the percentages 
found for ellipticals in previous works. 

3. The XLF of these sources is well fitted by a simple power law down to L x = 3 x 10 37 
erg s _1 (or to Lx = 2 x 10 37 erg s _1 after completeness correction), with a slope 
{(3 = 2.0 ±0.3) similar to that found previously for other ellipticals. The ratio between 
the total Lx due to LMXBs and the galactic Lb is lower than the average found for 
early-type galaxies, but still within the observed scatter. 

4. At the position of the galactic center a source (S2) is detected for the first time. Its 
spectral shape is quite hard (T = 1.49lg 13 ), without intrinsic absorption. S2 is however 
slightly extended, and a 0.3-8 keV upper limit of 2.8 x 10 38 erg s _1 is derived for a 
point-like nuclear source, one of the smallest established with Chandra. 

5. Three other sources, with a spectral shape typical of LMXBs, are detected in the 
central galactic region; they are as bright as the brightest sources in the XLF. Only 
one (S3) is consistent with being pointlike; the other two (SI and S4) are extended, 
and could be due to the superposition of few LMXBs and/or truly diffuse emission. 
The morphology of SI resembles a jet-like feature. 

6. Diffuse emission is detected out to a radius of ~ 30". A few independent lines of 
investigation, exploiting the spectral and imaging capabilities of Chandra, consistently 
indicate that this diffuse emission is due to unresolved LMXBs, and provide no evidence 
for hot gas, either close to the center or on a larger scale. The spectral analysis gives a 
(lcr) limit of L(0.3-8 keV)< 1.34 x 10 38 erg s _1 on any soft thermal component within 
a radius of 30". 

7. Other unresolved stellar sources (mostly M dwarfs and RS CVn systems) could con- 
tribute to the diffuse emission by < 10%. The general scaling of their collective Lx 
with the galactic Lb was also derived. 

8. Numerical simulations of the hot gas evolution for a galaxy model tailored on NGC 821 
show that the bulk of the gas is driven out in a wind for the whole galaxy lifetime, 
due to the heating provided by type la supernovae. While this may have been ex- 
pected based on simple energetic calculations, the simulations also show that the gas 
is accreting towards the center from within a very small inner region. This gas flow 
pattern is expected to be common in low Lb ellipticals, due to their cuspy central mass 
distribution. Gaseous accretion alone at the rates given by the simulations cannot have 
formed the nuclear MBH of NGC821. 
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In a companion paper (Pellegrini et al. 2007), we present further observational results 
for the central region of NGC 821, obtained with Spitzer and the VLA, and discuss various 
possibilities for the nature of the accretion process in this very low luminosity nucleus. 

We thank F. Ferraro for useful information related to Sect. I5.6[ the referee for useful 
comments and A. Kundu for kindly providing the list of globular clusters used in Sect. I3.1l and 
information on possible optical counterparts of the X-ray sources in Tab. [31 S.P. acknowledges 
partial financial support from the Italian Space Agency ASI (Agenzia Spaziale Italiana) 
through grant ASI-INAF 1/023/05/0. Partial support for this work was provided by the 
NASA Chandra Guest Observer grant GO5-6110X and by the Chandra X-ray Center NASA 
contract NAS8-39073. The data analysis was supported by the CXC CIAO software and 
CALDB. We have used NASA NED and ADS facilities, and have extracted archival data 
from the Hubble Space Telescope archive. 
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Table 1. NGC 821: main properties 



Type a 


B° a 


D b 


log L B 


Size a 


R e ° a e d 




M BH f 


1" 




(mag) 


(Mpc) 


(Lb, q ) 


(arcmin) 


(",kpc) (km s- 1 ) 


(cm -2 ) 


(1O 7 M ) 


(pc) 


E6 


11.72 


24.1 


10.27 


2.57x1.62 


43.9, 5.1 209 


6.2xl0 20 


8.5±3.5 


117 



a Type, Bj, and size from de Vaucouleurs et al. (1991; RC3). The size gives the major and minor 
axis of the D25 ellipse, that is the 25.0 B-mag/square arcsec isophote. The position angle is 25° 
(RC3). 

b Distance D from Tonry et al. (2001). 

c Effective radius R e in the R-band (from Soria et al. 2006b; see also Sect. 15.31) . 
d Effective stellar velocity dispersion (averaged over R e ) from Pinkney et al. (2003). 
c Galactic hydrogen column density (Dickey & Lockman 1990). 

f Gebhardt et al. (2003) report a value of 3.7^8 x 1O 7 M , later revised to the value given here 
(Richstone et al., |astro-ph/ 0403257) that is considered more reliable (Gebhardt, K. 2006, private 
communication) . 
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Table 2: Chandra ACIS-S observing log for NGC 821. 



Obs ID 


Date 


Exposure (ks) 


4408 


2002, Nov 26 


25.3 


4006 


2002, Dec 01 


13.7 


5692 


2004, Dec 04 


28.0 


6314 


2005, Jun 20 


40.1 


6310 


2005, Jun 21 


32.4 


6313 


2005, Jun 22 


50.1 


5691 


2005, Jun 23 


40.1 



Table 3. Properties of the X-ray point source population within D25 
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23.1 


_L1 n ^2 (S 


0.14 
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301.0 
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116.7 


8.10 


—0.44 


('—0 AG _0 ^Q 1 ) 


—0.27 
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(n 37 n ^n^ 


J020820.8+105909 b 


2 8 


20.8 


+10 59 09.9 


0.51 


32.01 


19.2 


8.6 


7.4 


0.52 


-0.70 


(-0.86 -0.37) 


0.48 


(0.14 0.94) 


0.09 


(-0.44 0.77) 


J020822.9+105959 


2 8 


22.9 


+10 59 59.1 


0.27 


32.30 


44.6 


8.7 


17.3 


1.20 


-0.62 


(-0.74 -0.45) 


-0.21 


(-0.47 0.03) 


0.55 


(0.37 0.92) 


J020823.3+105950 


2 8 


23.3 


+10 59 50.6 


0.30 


33.26 


43.3 


8.0 


16.8 


1.17 


-0.54 


(-0.63 -0.35) 


0.03 


(-0.16 0.20) 


0.18 


(0.04 0.39) 


J020819.1+110003 a 


2 8 


19.1 


+11 00 03.6 


0.47 


37.02 


11.4 


5.9 


4.4 


0.31 


0.10 


(-0.45 0.68) 


-0.01 


(-1.03 0.87) 


-0.35 


(-1.14 0.17) 


J020823.0+110008 


2 8 


23.0 


+11 00 08.9 


0.35 


39.19 


30.8 


7.7 


11.9 


0.83 


-0.31 


(-0.42 -0.02) 


-0.95 


(-2.23 -0.39) 


0.24 


(0.09 0.49) 


J020822.0+105904 b 


2 8 


22.0 


+10 59 04.6 


0.44 


39.42 


29.0 


7.5 


11.2 


0.78 


-1.00 


(-1.00 -0.73) 


-0.31 


(-0.67 -0.02) 


0.77 


(0.37 1.55) 


J020818.7+105919 b 


2 8 


18.7 


+10 59 19.3 


0.43 


41.70 


26.4 


7.4 


10.2 


0.71 


-0.56 


(-0.81 -0.29) 


-0.02 


(-0.34 0.27) 


0.52 


(0.13 1.21) 


J020819.9+105900 b 


2 8 


19.9 


+10 59 00.0 


0.54 


45.41 


17.1 


6.6 


6.6 


0.46 


-0.20 


(-0.61 0.21) 


-0.57 


(-1.67 -0.06) 


0.41 


(0.02 1.02) 


J020820.7+105855 a 


2 8 


20.7 


+10 58 55.4 


0.30 


46.64 


52.8 


9.2 


20.4 


1.42 


-0.46 


(-0.59 -0.32) 


-0.40 


(-0.68 -0.16) 


0.57 


(0.37 0.78) 


J020819.3+105900 d 


2 8 


19.3 


+10 59 00.2 


0.53 


49.23 


13.6 


6.1 


5.3 


0.37 


0.66 


(0.20 0.84) 


0.26 


(-0.65 1.34) 


-0.69 


(-1.81 -0.14) 


J020824. 1+110007 


2 8 


24.1 


+11 00 07.8 


0.31 


51.53 


43.6 


8.5 


16.9 


1.17 


-0.19 


(-0.35 -0.03) 


-0.56 


(-1.06 -0.22) 


0.26 


(0.13 0.47) 



Table 3 — Continued 



CXOU Name 




RA 


DEC 


A 


R 


cut 


err 


Fx 


Lx 


HR 


err 


C21 


err 


C32 


err 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


J020817.7+105907 a 


2 


8 17.7 


+10 59 07.7 


0.28 


60.28 


60.9 


9.6 


23.6 


1.64 


-0.72 


(-0.82 -0.60) 


-0.36 


(-0.54 -0.20) 


0.68 


(0.51 0.93) 


J020820.2+110041 d 


2 


8 20.2 


+11 00 41.1 


0.35 


61.17 


24.4 


7.2 


9.5 


0.66 


1.00 


(0.66 1.00) 


-0.30 


(-1.39 0.36) 


-0.49 


(-0.95 -0.15) 


J020822. 5+110044 


2 


8 22.5 


+11 00 44.1 


0.48 


65.73 


15.0 


6.4 


5.8 


0.41 


0.24 


(-0.29 0.72) 


-0.73 


(-2.06 -0.09) 


0.11 


(-0.27 0.57) 



Note. — Units of rigth ascension are hours, minutes and seconds, and units of declination are degrees, arcminutes and arcseconds. Units of A (the positional uncertainty) 
and R (the distance from the galactic center) are arcseconds. The positional uncertainty is based on the statistical error of the X-ray data alone (Sect.[3]|, while the absolute 
positional error is discussed in Sect. [4721 Units of Fx are 10 — 16 erg s _1 cm -2 and those of Lx are 10 38 erg s — 1 ; both quantities have been absorption-corrected for the 
Galactic column density and refer to the 0.3-8 keV band. See Sect.|3]for more details (e.g., on the definitions of the hardness ratio HR and the colors C21 and C23, and the 
calculation of their errors). 

a The position of this source is coincident with that of a candidate globular cluster found in the H ST /WFPC2 images, see Sect. 13.11 the offset between the optical and 
X-ray position is always well within A, except for the first of these 8 sources (J020821. 4+105946) that has an offset of 1.24A. 

b This is a field LMXB candidate (see Sect. 15m. 

C A specific spectral analysis has been done for this source (Sect. l3~Tl and Tab.|4]l. 

d This X-ray source may correspond to a faint background galaxy, from an inspection of the HST images discussed in Sect. |3~T1 (A. Kundu, private communication). 
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Table 4: Spectral analysis results for the brightest sources in Tab. [3j 



J020821. 5+105948 J020822.2+105922 J020823. 1+105952 



wabs(pow) : 

N„ (10 21 cm- 2 ) 2.2±25 0.5^:1 1.5j*i 

T 1.98j££ 1.21™ l-58±g:ii 

X 2 /dof 8.1/4 5.2/7 7.6/11 

Flux (10~ 15 erg cm" 2 s" 1 ) 4.5 9.7 10.9 

Luminosity (10 38 erg s _1 ) 4.1 6.9 8.9 



Note. — The Nh values are in addition to the Galactic one. Fluxes are observed values, luminosities are 
corrected for absorption, and are calculated for the 0.3-8 keV band. Errors give the 68% confidence interval 
for one interesting parameter. See also Sect. 13. ll 
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Table 5: Spectral analysis results for the central sources S1-S4. 



SI S2 S3 S4 



Net counts 178±14 246±16 324±18 82±9 
wabs(pow) : 

(10 21 cm' 2 ) 1.71°;* < 0.5 2.6+°| 3.01^ 

r 1 sn+ - 35 i /io+ - 14 o oq+o.23 <-> 90+0.52 

1 i.su_ 023 i.4y_ 013 /!.^_ 017 ^s_ 034 

rVdof 7.3/6 11.5/9 15.8/12 2.1/5 

Flux (10~ 15 erg cm- 2 s- 1 ) 5.8 8.6 9.1 2.3 

Luminosity (10 38 erg s" 1 ) 4.7 6.0 9.0 1.9 

wabs(apec) : 



(10 21 cm" 2 ) 


< 1.1 


< 0.5 




1 4+i-s 

L -^-1.3 


kT (keV) 


r 0+4.7 
°-°-1.8 




4.3+S:° 


3-4±g;| 


X 2 /dof 


7.0/6 


11.9/9 


18.4/12 


1.9/5 


Flux (10~ 15 erg cm -2 s -1 ) 


5.7 


8.1 


9.7 


2.4 


Luminosity (10 38 erg s' 1 ) 


4.3 


5.6 


7.6 


1.7 



Note. — The Njj values are in addition to the Galactic one. Fluxes are observed values, luminosities are 
corrected for absorption, and are calculated for the 0.3-8 keV band. Errors give the 68% confidence interval 
for one interesting parameter. 



Table 6. Spectral analysis of the diffuse emission (projected regions). 



Radius 


Net counts 


r 


F(0.3-8 keV) 
(10~ 15 erg cm -2 s _1 ) 


L(0.3-8 keV) 
(10 38 erg s" 1 ) 


kT 

(keV) 


F(0.3-8 keV) 
(10~ 15 erg cm -2 s _1 ) 


L(0.3-8 keV) 
(10 38 erg s" 1 ) 


X 2 /c!of 


5" 

10" 

30" 


166±13 
175±16 
581±48 


, 77 +0.17 

1 fi1+ - 23 
i - Di -0.21 
i 01 +0.13 

■■••^-o.ia 

i -°'-0.15 
1 74+0-14 

1 KO+0-09 
i - oy -0.17 


8.4±0.8 
8.6±1.2 
9.8±0.7 
10.0±1.4 
17.9±1.3 

i6.9il:J 


6.7±0.7 
6.7±1.0 
7.9 ± 0.6 
7.8 ± 1.1 
14.3 ± 1.10 
13 n+ 2 -40 


91+ - 16 

2^+ 08 
u - zo -0.06 

U - OM -0.18 


n Q7+13-6 

D RQ+0- 35 
u - DM -0.32 

1 06+ 042 

i.uo_ n 


44+ 16 2 
u - 44 -0.22 

n 7I ;+0.38 

n QR+0- 38 

u - yD -0.64 


8.6/9 
7.4/7 
25.0/19 
21.9/17 
78.2/64 
71.9/62 



Note. — Njj is fixed at the Galactic value of 6.2 X 10 cm" " 2 (Tab.[lJ; in the apec spectral model, the abundance is fixed at the solar value. 
Fluxes are observed, luminosities are intrinsic. 

Errors give the 68% confidence range for one interesting parameter. 
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Fig. 1 — Merged ACIS-S image of NGC821, with the D25 ellipse of Table Q] overplotted in 
green and the X-ray sources detected by wavdetect marked in red with circles of 3" radii 
(Sect. E]). 
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Fig. 2. — X-ray colors of point sources, left: with errorbars, and right: without errorbars 
(Sect. 13. ip . The grid indicates photon index values of 0, 1, 2, 3, 4 (from bottom to top) and 
N H values (from right to left) of 10 20 ,10 21 ,5 x 10 21 (dashed one), 10 22 cm . Red circles 
correspond to sources at R < 10", blue circles to sources between R = 10" and the ellipse 
D25. Filled circles are sources with net counts > 30, open circles with net counts < 30. 
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Fig. 3. — The Chandra ACIS-S spectra of the brightest point sources in Tab. [3] together 
with their best fit power law models (Tab. H]). From top to bottom: J020821. 5+105948, 
J020822.2+105922 and J020823.1+105952 (see also Sect. EH]) . 
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Fig. 4. — Cumulative XLF of point sources detected in the merged ACIS-S image within 
the D25 ellipse, with a central circle of 10" radius excluded (Sect. 13.21) . Black squares with 
error bar show the bias-corrected function; green circles are the uncorrected function; the 
red histogram gives the best fit single power-law; the dotted line shows the expected number 
of cosmic X-ray background sources from the logN-logS relation of ChaMP plus CDF data 
(Kim et al. 2007). 
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Fig. 5. — Co-added Chandra ACIS-S image of a central field of 25" x 25" in three bands: 
0.3-1 keV (top left), 1-2 keV (top right), and 2-4 keV (bottom left). North is up and east 
is to the left. The bar below the figure gives the correspondance between the color scale and 
the number of counts/pixel (the size of 1 px is Of! 5). See Sect. H]for more details. 



Fig. 6. — Co-added smoothed "true-color" image of a central 55" x 55" field (see details 
in Sect. 0J. The red, green and blue colors correspond to the 0.3-1 keV, 1-2 keV and 2-4 
keV bands respectively. North is up and east is to the left. The diffuse emission follows the 
optical shape of the galaxy (see Tab. [fl Fig. [I] and Fig. [8] below). 
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Fig. 7.— The ACIS-S image in the 0.3-6 keV band of the central 40"x25" region of NGC 821, 
with the sources and their extent detected by the CI AO task wavdetect (Sect. 14. ip . North is 
up and east is to the left. Here 1 px=0'/25. The semi-major and semi-minor axes for S1-S4 
are respectively 2"5 x 1"5, 1"7 x 1"5, V!2 x 1"2 and V!2 x l'/O. 
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Fig. 8.— HST WFPC2 image of NGC821, with the positions of the Chandra ACIS-S 
sources detected by wavdetect marked with red circles. The six yellow circles show the best 
optical/X-ray coincidences falling within the D25 ellipse (not including the nuclear source) 
that were found in Sect. 13.11 and were also used for the astrometric analysis of Sect. 14.21 
HST isophotes for the central galactic region are also overlaid. 
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Fig. 9. — ACIS-S image of NGC821 with green circles of 0'!7 radii marking the positions 
of the CIAO wavdetect sources (lpx=0'/5), and the yellow circles corresponding to GCs (as 
in Fig. [8]). The galactic isophotes are also overlaid. In the recalibrated WFPC2 image, the 
galactic center is located at RA=02 h 08 m 21113, Dec=+10° 59' 41'.'8 (see Sect.S^D and falls 
within the circle of the Chandra source S2. 




Fig. 10. — The Chandra ACIS-S spectra of the extended sources S1-S4 together with their 
best fit power law models (Tab. EJ). From left to right the panels refer to SI and S2 (top), 
S3 and S4 (bottom). See Sect. 14.31 for more details. 
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Fig. 11. — Radial profile of the (0.3-6 keV) emission, after removal of sources detected by 
wavdetect (Sect. loTTj) . The flattening at radii > 25" is due to field background. Vertical bars 
give the ±1<t uncertainty. 
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Fig. 12. — Radial profiles of: detected point sources (in blue), in units of number /square 
arcmin; background-subtracted diffuse emission over 0.3-6 keV (in red), in units of 
counts/pixel; galactic R-band emission (dashed lines, arbitrarily normalized). See Sect. 15.31 
for details. 
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Fig. 13. — Radial distribution of the diffuse 0.3-1.5 keV emission (in black) that is implied 
by the observed 1.5-6 keV emission, when assuming for it the spectral model described in 
Sect. 15.41 In red the observed profile in the same 0.3-1.5 keV band is also shown. The 
flattening at large radii is due to field background. Vertical bars give the ±lcr uncertainty. 



